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Abstract. In the very early Universe the matter may be described by the free radiation, 
that is by the set of massless fields with negligible interactions between them. Then the dom- 
inating quantum effect is the trace anomaly which comes from the renormalization of the 
^ \ conformal invariant part of the vacuum action. The anomaly-induced effective action can be 
found with accuracy to an arbitrary conformal functional which vanishes for the special case 
of the conformally flat metric. This gives the solid basis for the study of the conformally-flat 
cosmological solutions, first of which was discovered by Mamaev and Mostepanenko and 
by Starobinski in 1980. Treating the anomaly-induced action as quantum correction to the 
Einstein-Hilbert term we explore the possibility to have inflationary solutions, investigate 
their dependence on the initial data and discuss the restrictions in considering the density 
perturbations. The shape of inflationary solutions strongly depends on the underlying gauge 
model of the elementary particles physics. Two special cases are considered: Minimal Stan- 
dard Model and the matter sector of N = 8, D = 4 supergravity. It turns out that inflation 
is almost inevitable consequence of the great difference between Planck mass and the mass 
of the heaviest massive particle. 
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1 Introduction. 

The standard cosmological model JTJ was constructed soon after the creation of General 
Relativity and it describes the present stage of the Universe pretty well (see, for example, 
At the same time this model fails to explain the early period of the evolution, and 
hence some new theoretical idea is required. It is commonly accepted that problems of initial 
singularity, horizon, flatness and monopoles may be solved if one supposes that the early 
Universe passed through the period of a very fast (more than the speed of light) inflation when 
its size increased in more than 60 orders. Moreover, the origin of cosmological perturbations 
has an elegant explanation within such scenario. On the other hand, the source of inflation 
still remains unclear. In the first works devoted to the exponential expansion of the Universe 
the origin of it was considered in the quantum effects of vacuum ^ |J (see also ||), but the 
real explosion of interest to it started after the paper M where inflation was induced as a 
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consequence of the Spontaneous Symmetry Breaking in the Standard Model of electroweak 
interactions. The further development has shown that the last cosmological scenario was 
not perfect and many other models, theories and approaches were considered. Finally it 
was realized that all the requirements to the inflationary model can be satisfied only if one 
introduces some new scalar field called inflaton, which has a potential of a very special form 
(see, for example, [0, [| for a comprehensive review). 

Here we intend to go backward to the historic origin of inflation []3], |J and derive it from 
the quantum effects of pure radiation on curved background. It will be shown that the 
inflationary solution doesn't require any fine-tuning or special choice of the parameters of 
the underlying conformal invariant matter fields but is absolutely natural and exists for a 
wide set of initial data. The exit from the inflationary phase may be also naturally explained 
by the nonstability of the exponential solution with respect to the density perturbations that 
can be relevant at the end of inflationary epoch. In this letter our consideration is restricted 
by the conformally flat metrics, but this does not mean that the similar results can not be 
achieved for the k = +1 or k = — 1 cases. 

The paper is organized as follows. In the next section we consider the model for the 
very Early Universe filled by the free radiation. This radiation is composed by the set of 
free, massless, conformally invariant matter fields on an arbitrary curved background. We 
review the derivation of the quantum correction to the classical action of gravity [[U], |TT 



This correction results from the trace anomaly |T2| and it contains an arbitrary conformal 
invariant functional. However, for the conformally flat background metric this functional is 
irrelevant and hence in the framework of the above model the solution of [|K], |TTJ is exact. In 
section 3 the particular inflationary solution is presented and its nonstability is established. 
Section 4 contains the results of numerical investigation and the physical analysis of the 
inflation. In the last section we draw some conclusions and outline the prospects for further 
studies. 



2 Free fields and quantum corrections 

The action of the ideal liquid filling the Universe can be written as 

S M = J d 4 x^j | -- (p + p) u»u u g^ + ^p + -pj . (1) 

It is easy to see that for the conformally flat cosmological metric the variation of this action 
with respect to the conformal factor of the metric 

8Sm _ , Q n 3 
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vanishes for p = |, and therefore = , so the radiation should be described by the 
conformal invariant field actions. 

In the framework of the asymptotically free gauge theories the interactions between 
matter fields are weakened in the high energy limit. Therefore we start from the action of 
free massless conformal invariant fields: No scalars (spin-0), Nu^ spinors (Dirac, spin- 1/2) 
and N\ abelian vectors (spin-1). All iV's indicate a number of fields (not multiplets) in 
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curved space-time, taking conformal coupling with scalar curvature for the scalars. 



S„ 



matter 
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For the free fields on curved background the only possible divergences are the one-loop 
vacuum ones. Therefore the renormalizability requires the action of vacuum to be introduced 
in the form (see [13 and [[14]] for the introduction) 

Svacuum = J d*X^/=g~ {hC 2 + l 2 E + l 3 nR} (3) 



where h,2,3 are some new (with respect to the flat space-time) parameters, C 2 is square of the 
Weyl tensor and E is the integrand of the Gauss-Bonnet topological (in d — 4) invariant. The 
vacuum action (|3|) includes only conformal invariant and surface terms. For the convenience 
of the reader we reproduce known formulas 
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One has to notice that the introduction of conformal non-invariant terms like / ^/—gR or 
/ y/—gR 2 is possible but not necessary for the renormalization of the free conformal invariant 
theories. The renormalization of vacuum leads to the conformal (trace) anomaly |12] which 
enables one to find the nonconformal part of the effective action of vacuum |TT|. For the 
sake of completeness and also to fix notations and signs we reproduce the main steps of this 
derivation. 

The form of the one-loop divergences is well-known (see [0, |T4|j). 
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The contribution of the Weyl spinor is half of the Dirac one. We remark that the last two 
(topological and surface) terms in (f|) are very important because they contribute to the 
trace anomaly and thus affect the dynamical equations for the effective action of gravity. 
The one-loop effective action has the form 



s + t, r = r^ + A5 



(6) 
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J J ± 05 



3 



where I^ 1 -* is naive quantum correction to the classical action S and AS is a counterterm. 
The anomalous trace of the Energy- Momentum Tensor is |12| 



T =< >-- 



wC 2 + bE + cUR 



(7) 



with the same coefficients a,b,c. The Eq. (0), in turn, gives rise to the equation for the 
finite part of the 1-loop correction to the effective action 
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The solution of this equation is straightforward [10], |TT| and one gets, in terms of g^ u 
9fiv ■ exp[2a], the following quantum correction to the vacuum action []. 

f = S c [g^\ + f d 4 xy/=g~ {waC 2 + ba{E- -OR) + 26a Aa- 



-^(c+lb)[R-Q(Va) 2 -6(na)] 2 )}. (9) 

Here the fiducial metric g^ has fixed determinant, A is conformal invariant self-adjoint 
operator 

A = a 2 + 2BTVJ7 V - \ru + ~ (WR) 

and S c [g^ v ] = S c \g^ is some unknown conformal-invariant functional. In general, exact 
calculation of this functional is impossible. However, if we are interested in the conformally 
flat cosmological solutions, this functional is (fortunately) of no importance for us. 



3 Particular inflationary solution 

Our purpose is to implement the quantum correction ([|) into the classical gravity-matter 
system. For the sake of simplicity here we consider only the conformally flat metrics, for 
which S c [g^ y ] is irrelevant. The solution (|9|) should be added to the classical action 5^0^^ + 
Smatter- However, since we are interested in the conformally flat matter, it decouples from the 
conformally flat metric 0. Hence the only nontrivial contribution to the dynamical equation 
for a comes from the vacuum part. The conformal and surface terms in the vacuum action (|3]) 
do not contribute to the dynamical equations for the conformal factor. On the other hand, 
in order to have correspondence with classical gravity one has to add the Einstein-Hilbert 
term to the vacuum action. Then we meet the total action of the form 

S t = g / d^x^f—g R + T[g] + conformal invariant terms. (10) 

5 Recently there were indications Jl^, |l6| that the solution (0) does not assume the test based on the 
calculation of the three-point function, and hence another solution for Y should be found. We emphasize that 
any solution for the effective action differs from (^|) by the conformal invariant functional and this difference 
is not relevant for the cosmological solutions which will be explored here. 

6 This fact can be seen, in particular, from the identity T£ = which holds for the Energy-Momentum 
Tensor of the radiation. 
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and, denoting conformal time as 77 and taking conventional a = In a(j]) we arrive at 
/ d^x^J—g R = 6 / d 4 x a' 2 . Here the prime stands for the derivative with respect to 77 . 

Let us remind that kT 2 = Mp t . Taking variational derivative of the total action with 
respect to 0,(77) one arrives at the equation of motion for 0.(77): 



a"" a' a"' a" z / 2b\ a" o' z 46 o /4 2 




-4— __3-=- + 2 [3--]—- + -— --M 2 Pl a"a = 0. (11) 



c 

It is convenient to rewrite ([11]) in terms of physical time t where a(rj)dr] = dt so that 
d/drj = ad/dt etc, and denote the derivatives with respect to t by points. 

2 • - ( 46\ .2 .. ..2 2.1/r,, / •■ -2\ 
a a+3aaa — I 5 H j a a + aa [a a + aa J =0. (12) 

The last equation contains two dimensionless constants b, c which depend on the particle 
content of the matter fields and may be fixed by choosing particular gauge model. Also there 
is a dimensional constant Mpi , which defines the scale. The Eq. ( |T2"D looks quite complicated 
and it is difficult to solve it exactly. Here we explore a particular exponential solution and in 
the next section perform the numerical study and physical analysis. Below we present all the 
numerical results for the Minimal Standard Model and and for the compactified to D = 4 
supergravity (M-theory). Indeed, since gravity is treated as purely classical background, in 
the last case one has to disregard all spin-2 and spin-3/2 degrees of freedom and keep only 
contributions of the spin-(0, 1/2, 1) sectors. 

Let us look for the exponential solution of the form a(t) = aa(t) — A ■ exp(Xt) where 
A, A are some constants. Substituting this expression into Eq. fll2]) one can easily see 



that it gives solution for an arbitrary A and A = ±77=5- Hence the crucial point is that 
b in (pp is negative for any particle content of the original theory (0). There are always 
two solutions with opposite signs for A , and the positive one describes the exponentially 
expanding Universe. This exponential solution is exactly the one which has been obtained 
in [|], U by the use of the renormalized Energy- Momentum Tensor []T7| (see also 0]). Some 
comment is in order. The equation for a(t) which was used in [|J looks different from (|P2|). 
This is because the in M the (0, 0)-component of the equation 



R^u - g Rg^u = 8tt < > 



has been used while we use the trace of it. Indeed both equations contain the same informa- 
tion since the metric has only one degree of freedom. On the other hand, in our framework 
it is clear that the introduction of the conformal invariant matter fields (or the action for 
ideal liquid ([!]) with p — |) doesn't change the equation, so the above inflationary solution 
is valid for any (conformal) matter content]]. An important consequence is that the density 
perturbations which do not violate the constraint p = | can not destroy the exponential 
behavior of a(t) . On the other hand, any perturbations which do violate this constraint 
imply the appearance of the massive particles or cosmological constant. In both cases one 
can not consider these perturbations in the framework of the above model based on the 
conformal invariant fields (0) and conformal anomaly (0). If the initial theory contains some 
massive fields, than the trace < > has non-anomalous (both classical and quantum) 

7 In this point our analysis differs from the one of [^] . 
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contributions and one has to take them into account. Suppose the initial matter fields have 
masses of the order of 100 GeV that is typical for the Standard Model. Then this matter 
can be successfully considered as (fj) at the energies a few orders below the Planck scale 
10 19 GeV and then the Universe is inflating until the typical energy becomes comparable 
with the masses. After that the density perturbations violating the constraint P = § become 
relevant and they provide the successful exit from the inflationary phase and also basis for 



the formation of the galaxies ||18|| . Of course this scheme implies that below a Planck mass 
there is a sufficient gap on the mass scale, and the heaviest massive particle below M p i is 
many orders lighter. If this condition is satisfied, one can regard all particles as massless, so 
that the constraint p = p/3 takes place. Then the most important questions concerning the 
exponential solution above are: 



1. whether the inflationary solution is stable under the metric perturbations which are 
consistent with p = |? 

2. whether the period of inflation is sufficient to provide the necessary rate of expansion? 



3. whether the inflation occurs only for the very special initial data or it can be achieved 
in a more general situation? Indeed the above solution a(t) = A ■ exp(Xt) implies that 
a (™)(0) = A n a(0), . Hence this solution is only valid for some 1-dimensional line in the 
4-dimensional space of the initial data. 

In order to understand better the sense of the Planck mass in the above expressions, one 
can introduce the dimensionless time r = Mp\ ■ t . The equation (|l|) remains almost the 
same, the only change is Mp\ — > 1. Then, in terms of the Planck units we have the following 
inflationary solution 

a (t) = Ae +T/V ~ b (13) 

The Eq. (|T3|) can be used to calculate the necessary duration of inflation. Suppose we 
want the Universe to expand in 10™ times, starting from some moment To. According to 
([131) , the total rate of inflation is 



a ° (T0 + Ar) expiW^Arl, AT, = ^ + 11.^ + 62-^, (14) 



a (r ) " [ V N t 

and thus we arrive at 



In 10 ., 
Ar = — — t/— • n. 




For the Minimal Standard Model JV = 8, JVi = 48, N\ = 12 and we arrive at Ar « 0.345-n . 
Taking n = 60 one meets the necessary period of inflation Ar rs 20.7 1 pi where tp\ is the 
Planck time, which serves as a time unit for r. Hence the necessary period of inflation is 
just one order higher than the Planck time. If one increases the particle content of the 
theory Ar becomes longer. In particular, this happens when we use SU(5) GUT or the 
supersymmetric models. One can remind that the contributions of different fields to the 
coefficient b in (0) have the same signs. In general case the expected time of inflation can 
be evaluated as one-two orders longer than the Planck time, depending on the choice of the 
model. 
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As it was already mentioned above, it is important to investigate the metric perturbations 
around the solution (|i~3l). However this meets serious difficulty. While taking an arbitrary 
perturbations one can not avoid variations of the fiducial metric g^ v . In this case one has 
to study this perturbations on the basis of the full action ( |10"D and the conformal invariant 
terms may be as important as the anomalous correction which we deal with. In principle 
one can mimic the conformal action by the expressions described in [[TjJ but we prefer to 




postpone this consideration for some future work. Here we consider only the perturbations 
of the conformal factor. It proves useful to rewrite the equation ([12]) for the cr(r) = In a(r). 

= (15) 

and only then perform the perturbations of er(r) = cro( r ) +x(t) near the inflationary solution 
o"o(t) = ha. A + -^^t ■ For the sake of simplicity below we put hi A — 1. The equation for 
the perturbations have the form 

7 ... 2 / b\ .. 4 . 
x H — -= x — - 6 x -== x = (16) 



b 



(all derivatives with respect to r) and can be easily explored. The stability of the inflationary 
solution under the perturbations of conformal factor strongly depends on the particle content 
No, Nu2, N± of the initial action (^j. We remark that the stability depends only on the ratio 
b/c. With the proportionally increasing number of fields iVo,i/2,i the necessary period of 



inflation becomes longer but the (non) stability of solution (13) doesn't change. 

Substituting x = e r ' T one meets three different nonzero roots r lv , 3 . We remark that 
the existence of the zero root manifests the arbitrariness of one of the coordinates in the 
4-dimensional space of initial data (choice of In A). This can be also seen if one makes 
perturbations in ( |12"D but in the last case a bit more sophisticated consideration is required. 
The values of r lj2 ,3 depend on the particle content AT 0)1 / 2 ,i ■ Taking the values typical for the 
Minimal Standard Model we arrive at 

n /2 « -9.99 ± 24.71 i , r 3 « -26.66 . (17) 

Thus in the case of Minimal Standard Model all four roots have negative real parts. This 
means that the exponential inflation is stable under the small perturbations of the conformal 
factor of the metric. The possible sources of instability remain in the perturbations of the 
other degrees of freedom or in the density perturbations discusses above. It is interesting 
that the stability analysis performed (in a different way) in Q gave different result. The 
source of the discrepancy is that in |4[] the negative sign of c in ([7]) has been chosen. This 
sign, indeed, corresponds to the contributions of vector fields into @, while spin-0 and spin- 
1/2 fields contribute with positive sign. Here we presented the results for the field content 
of the Standard Model, and will confirm our conclusion using numerical methods. 

Let us consider, as a second example, the M-Theory, which is nowadays regarded as a 
possible candidate for the unified model of all interactions. In the framework of M-theory one 
possible realization is maximal N = 1, D — 11 supergravity. In this case its particle content 
is unique and the gauge group is SO (8). Then the spin-0, 1/2, 1 sector of the compactified 
to D = 4 theory has 70 real scalars, 28 vectors and 56 Maiorana spinors. Taking iV(o,i/2,i) = 
(70, 28, 28) one meets the roots 

r i /2 Rs -54.63 ±71.39 i, r 3 sa +72.63 , (18) 
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indicating to the instability of the exponential inflation. The question is whether this means 
instability of inflation at all? As we shall see in the next section, the answer is exactly 
opposite. 

4 Numerical study and solutions with bounce 

The fourth order equation ([12]) depends on two parameters b and c, the first being negative 
for any matter content, and c being positive or negative depending on the particle multiplet. 
Numerical integration of this equation depends on the value and sign of b and c and on the 
initial conditions. Among the many possible scenarios, we can identify three main ones: 

1. Singular expanding Universe, with oscillations between accelerating and desaccelerat- 
ing phases; 

2. Inflationary expanding Universe; 

3. Non singular Universes, presenting a contraction initial phase, followed by an expand- 
ing phase. 

The exact inflationary solution ([l3| ) represents a particular choice of initial conditions. 
The perturbative study performed before indicates that this solution is stable for the particle 
content of the Standard Model (see Figures 1,2). 

Since it is inconvenient to make a numerical study (and especially graphical presentation) 
of the Universe suffering from 60-orders inflation, we have explored the behavior of logarith- 
mic variable cr(r) rather then the a(r) itself. It turns out that the inflation occurs for a 
very wide set of initial conditions, even when these conditions are taken far away from the 
exact exponential solution fllTf ). In some cases, depending on the multiplet composition, the 
stable solutions correspond to some kind of "hiperinflation" and expansion perform much 
faster than exponentially. 

A numerical study of the initial conditions around those representing the solution ( |T3"D 
reveals that this inflationary solution is "robust": there exist two main scenarios: 

1. An eternal inflationary expansion; 

2. A non singular Universe, presenting a bounce for the scale factor. 

In the last case, the Universe exhibits an accelerating contraction followed by an accelerating 
expansion; just in some special cases there is a short non inflationary period of contraction 
(expansion). These numerical results confirm the analytical one, and extend the possible 
scenarios including Universe exhibiting bounce. Only these two scenarios were found for the 
Standard Model (Figures 1,2 give typical examples of the plots). 

We remark that the cosmological scenario exhibiting bounces may be encounter in many 
non linear curvature models and in some special scalar-metric gravity theories. The main 
distinguishing feature of the scenarios obtained here is its almost universal inflationary be- 
havior, while in the other theories there was just a short period of inflation near the bounce. 
The continuosly expanding and bouncing solutions obtained here have the number of e— fold 
sufficient to solve the problems that can be cope with an inflationary phase. 

As a second example we consider the multiplet of the M-Theory described in the previous 
section. In this case, if one performs a slight change of the initial conditions around those 
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typical for the exact inflationary solution, there can be an extremely fast expansion, with 
singularity just a fraction of (Planck) time unite after the initial moment. Some "typical" 
plots are presented at the Figures 3,4. One can call the behavior shown at Fig. 3 as 
" hyperinflation" , while the one on Fig. 4 corresponds to some drammatic contraction of the 
Universe. Making the time inverse, however, we meet inflation also in this case. Solutions 
with the bounce are also possible. 



5 Conclusions and discussions 

We have considered the quantum effects of free radiation on the background of conformally 
flat metric and found that this leads, in a very natural way, to the inflation which must 
continue only a very short time to resolve the problems which inflation is supposed to resolve. 

It is important to notice that the existence of exponential solution and the velocity of 
expansion depends only on J E-type topological counterterm and absolutely doesn't depend 
on the J OR and J C 2 -type counterterms (|4]). This may be important, since in the same 
manner as we have added / R to the vacuum action of the conformal theory, one could also 
add local (classical) / i? 2 -term with an arbitrary coefficient which can mixture with the same 
term in (|9]). That is why the / i? 2 -based inflation is not safe in general. But in the case of 
flT3"D everything depends only on the nonlocal term in the effective action (0) which does not 
suffer from this ambiguity. The only physical requirement which is important for inflation 
is a gap on the mass scale between Planck mass M p i and heaviest massive particle M max . In 
the course of inflation the typical energy scale decreases and, when is becomes comparable 
to M max , the mass of the particles can be seen, the violations of p — | constraint show up 
and the inflation quits due to the perturbations. 

The remaining open question is indeed the stability of inflation with respect to arbitrary 
perturbations of the metric. As it was mentioned above, this question can not be addressed 
within the exact solution (|9]) because it is necessary to use some approximation (or imitation) 
for the unknown conformal invariant functional S c . Similar situation takes place for the 
unconstrained density perturbations which can occur at small energies (or, in other terms, 
at later period of the evolution of the Universe). In this case one has to derive, in some 
reasonable approximation, the effective action of the massive fields and, as a best option, 
also take into account their interactions. In case of the weakly interacting conformal field one 
has to count the effect of back reaction of vacuum to the matter fields which produces a slow 



running from the conformal fixed point p0[ . In this case one can explore the possibility to 
have inflationary solutions for the induced gravity. We are going to reconsider these problems 
elsewhere. The present letter contains a brief report on the possibility to have consistent 
and natural inflation from the vacuum effects and some review of its main features. 
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FIGURE CAPTIONS 




Figure 1: Typical inflationary behavior for the Standard Model case. Initial data: 
a(0) = 1, (j(0) = 2m/ v /z 6, <r(0) = 0, o=(0) = 0. 
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Figure 2: Another typical inflationary behavior for the Standard Model case. Initial data: 
a(Q) = 1, (j(0) = 0.1 • m/y/-b, <r(0) = 1, <r(0) = 0. 
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Figure 3: Typical inflationary (hyperinflationary) behavior for the M-theory case. Initial 
data: tr(0) = 1, &(0) = m/V^b, cr(0) = 10, o=(0) = 0. 




Figure 4: Non-inflationary behavior for the M-theory case. This kind of scenario with 
"hyperinflation" is possible only for large negative <j(0). Initial data corresponding to this 
plot: tr(O) = 1, <j(0) = m/y/-b, <r(0) = 0, o=(0) = -400. The change of the signs in o=(0) 
and (j(0) is equivalent to the inverse of time and leads to "hyperinflation". 
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